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Introduction 


Marsaglia  [1]  has  given  a  simple  method  for  generating  exponen¬ 
tial  random  numbers  on  a  digital  computer.  We  present  a  similar 
method  for  generating  random  numbers  with  the  chi  distribution.  Such 
random  numbers  may  be  used  to  generate  normal  random  numbers. 

I.  The  chi  distribution  (of  rank  two)  F  is 

F(a)  =  0,  a  <  O, 

F(a)  =  1  -  e~®  O  <  a. 


Let  X  be  a  random  variable  with  the  distribution  F.  Let 
G^(a)  =  Prob(x  <  a|x  <  c)  where  0  <  c.  Then  for  0  <  a  <  c, 

G^(a)  =  (1  -  e"®^/2)/(l  - 


2/  2.k 

=  1  -  E  q  (1  -  a  /c  )  , 
k=l 


where 

2 

q.  =  (c^/2)Vi:k:(e‘'  -  1)]. 

(  2  2./ 

Let  H(a)=l-e  'c<a, 

=0,  a  <  c. 

2  2 

Then  F(a)  ^  (1  -  e"°  /^)G^(a)  +  e"'^  /^H^(a). 


Thus  a  random  number  with  the  distidbution  F  may  be  generated 

6is  follows.  Generate  a  uniform  random  number  u,  i.e.,  a  random 

-0^/2 

number  uniformly  distributed  on  (0,1).  If  u  <  (1  -  e  ), 


2 


generate  a  random  number  wita  the  distribution  G^;  otherwise 
generate  one  with  distribution 

A  random  number  y  with  the  distribution  may  be  generated 

by  setting  y  =  t/a.r  +  c^,  where  r  is  a  random  number  with  the 
exponential  distribution. 

A  random  number  x  with  the  distribution  can  be  generated 

by  setting 

X  =  c  .  minCmaxCuj^.u^), . . .  ,max(u2^_j^,U2^)D, 
where  the  u^  are  independent  uniform  random  numbers,  6uid  z  is  a 
random  integer  taking  on  the  value  k  with  probability  q  .  This 

iC 

fact  is  easily  verified  by  noting  that  the  distribution  of  x  is 
just  the  series  (1). 

For  a  binary  computer  the  best  choice  for  c  is  c  =  2,  On  the 
IBM  7090  computer  the  average  time  to  generate  a  chi  random  number  x 
by  this  method  is  112  cycles,  (A  cycle  is  2.14  microseconds  on  this 
computer).  This  assumes  that  the  exponential  rsindom  numbers  are 
generated  by  the  method  given  in  £1]. 

If  X  is  generated  by  setting  x  =  the  average  time  is 

165  cycles, 

II.  To  generate  normal  random  numbers  we  make  use  of  the  following 
well-known  fact.  Let  (a,P)  be  the  rectangular  coordinates  of  a 
random  point  uniformly  distributed  on  the  unit  circle.  Then  if  x 
is  a  chi  random  number  y  =  ox  and  z  =  Px  are  independent  standard 
normal  rauidom  numbers. 
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The  following  methods  for  generating  such  a  pair  (a,^)  are 
well  known. 

Method  1.  Test  independent  pairs  of  uniform  numbers  (u,v)  until  a 

2  2  y  2  2 

pair  is  found  which  satisfies  u  +  v  <1.  Then  set  a  =  u/vu  +  v 
and  p  =  +  v^. 

Method  2.  Test  independent  pairs  (u,v)  until  a  pair  is  found  which 

2  2  2  2 
satisfies  u  +  v  <1.  Then  set  o  =  2uv/(u  +  v  )  and 

-  .2  2...  2  2. 

P  =  (v  -  u  )/(u  +  V  ). 

To  generate  normal  random  numbers  we  can  use  the  following  pro¬ 
cedure.  Generate  a  chi  random  number  x.  If  x  <  c,  use  method  2  to 

generate  (a,p).  If  c  <  x  use  method  1  to  generate  (a,p).  Note 

/2 

that  we  can  decide  if  c  <  x  before  we  set  x  =  vc  ♦  2r.  Therefore 
this  square  root  operation  can  be  combined  with  that  used  to  generate 
(a,P).  In  effect  when  c  <  x,  we  compute  a  pair  of  normeG.  random 
numbers  y  and  z  by 

y  =  xa  =  u[(2r  +  c^)/(u^  +  v^)]^ 

and 

z  =  xp  =  vC(2r  +  c^)/(u^  + 

This  procedure  takes  156  cycles  to  generate  one  normal  random 


number  on  the  7090. 
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